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Abstract  (continued) 

The  overall  goal  of  the  research  program  is  to  develop  an  atomistic  approach  to  gain  an 
understanding  of  the  mechanisms  of  growth  processes  and  to  contribute  to  the  development  of 
meta^'trietal  and  metal^semiconductor  heterostructures.  The  Phase  I  research  involved  (i)  reliable 
modelling  of  underlying  atomic  interactions  within  the  atomic  constituents  of  the  substrate, 
interface,  and  adlayer;  (ii)  static  and  dynamical  studies  of  interfacial  energetics  and  kinetics.  A 
survey  of  available  modelling  schemes  has  been  made  and  a  strategy  for  our  own  future  modelling 
efforts  is  identified.  y  COf>f>pr ,  fyuyi-  v\  <£l - 

v 

An  analysis  of  the  various  modelling  schemes  for  adatom  diffusion  studies,  based  on  the 
application  of  dynamical  calculations,  has  shown  that  (a)  all  the  models  exhibit  the  experimentally 
observed  directional  anisotropy  for  the  self  diffusion  on  Ni  (110)  surface,  (b)  estimates  of  the 
classical  diffusion  barrier  for  Ni/Ni  (100),  Ni/Ni  (110),  and  Ni/Ni(lll)  indicates  that  diffusion 
will  occur  more  readily  on  (1 1 1)  surface  orientation,  and  (b)  at  high  temperatures,  the  models  agree 
with  each  other  but  as  temperature  decreases  the  predicted  diffusion  data  fall  within  a  scatterband. 
The  latter  observation  points  out  the  necessity  for  precise  measurements  of  self-diffusion  of  Ni  on 
Ni. 

Atomistic  static  energetics  of  a  Ni  monolayer  on  a  Cu  substrate,  based  on  both  a  simple 
pairwise  additive  potential  and  a  more  realistic  many-body  potential  model,  indicate  that  the 
characterization  of  the  details  of  axial  commensuration  in  terms  of  energy  minima  and  maxima  is 
found  to  depend  on  the  atomistic  model  such  as  rigid  vs  nonrigid  (nonrigid  approximation 
corresponds  to  the  inclusion  of  adatom-adatom  interactions)  for  the  adlayer.  A  successful 
determination  of  the  separation  of  the  total  adatom-substrate  interaction  energy  in  the  realistic  model 
led  to  energy  charactersitics,  within  the  rigid  lattice  approximation,  that  are  similar  to  those 
obtained  with  the  simple  pairwise  potentials.  This  implies  that  for  rigid  lattices  the  energetics  are 
only  determined  by  the  relative  atomic  sizes  of  the  substrate  and  monolayer. 

The  strain  energy  for  Cu/Ni  is  larger  than  that  of  Ni/Cu.  This  indicates  that  a  large  number 
of  layers  will  grow  pseudomorphically  in  the  case  of  negative  misfit  interface,  Ni/Cu.  Studies 
involving  the  adatom  diffusion  and  the  strain  energy  dependence  on  the  substrate  orientation  may 
provide  a  clue  to  the  choice  of  substrate  orientation  for  obtaining  better  quality  films  and  for 
understanding  the  competition  between  energetics  and  kinetics  at  various  temperatures. 
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Introduction 

An  atomistic  understanding  of  the  metal-metal  and  metal-semiconductor  interfaces  is  highly 
desirable,  both  for  fundamental  physical  reasons  and  due  to  increasingly  sophisticated  technical 
applications.  Both  stable  and  metastable  phases  through  the  growth  of  epitaxial  thin  films  on 
suitable  substrates  have  been  observed  recently 1  and  are  leading  to  the  fabrication  of  superlattices. 
Although  many  of  the  growth  parameters  are  controllable  in  the  modem  molecular  beam  epitaxy 
(MBE)  and  metalorganic  chemical-vapor  deposition  (MOCVD)  techniques,2*3  the  precise  roles  of 
not  only  the  various  growth  parameters  but  also  of  the  energetics  and  kinetics  in  affecting  the 
diverse  phenomena  occurring  at  growth  fronts  are  not  known.  Furthermore,  sometimes  the 
precise  interpretation  of  the  nature  of  growth  becomes  difficult.  For  instance,  there  is  some 
controversy  in  the  recent  literature4  regarding  epitaxial  fee  Fe  films  on  Cu  (100).  A  microscopic 
approach  via  dynamical  simulations  for  such  issues  would  play  a  dominant  role  in  clarifying  this 
controversy  of  the  heterostructure  in  terms  of  the  energetics,  structure  and  kinetics  of  the  atoms 
involved. 

Although  there  have  been  a  few  atomistic  energetic  studies,5  these  studies  have  been  based 
on  the  Lennard-Jones  form  for  atomic  interactions  which  is  inadequate  for  metallic  and 
semiconductor  modelling.  These  preliminary  level  calculations,  however,  provided  a  clue  on  the 
strains  which  could  be  developed  in  the  deposited  layer  due  to  the  lattice  mismatch.  As  realistic 
interaction  potentials  are  being  produced  for  metals6  and  semiconductors,7  it  is  timely  that 
quantitative  information  concerning  growth  and  interface  be  obtained  by  adapting  these  potentials 
in  computer  simulations.  The  recent  intense  studies  on  homoepitaxial  growth  of  Si8  have 
addressed  a  number  of  interesting  issues:  (a)  the  growth  dependence  on  the  surface  orientation,  (b) 
growth  characteristics  as  a  function  of  temperature,  (c)  significance  of  annealing  after  deposition 
(d)  the  role  of  atomic  collisions,  (e)  atom  versus  cluster  deposition  and  (f)  growth  mode.  These 
studies  are  not  only  model-based  but  also  do  not  take  into  account  the  electronic  structural  features. 

First  principles  energy-band  calculations9  are  being  performed  on  these  new  materials  but 
they  do  not  include  finite  temperature  effects.  An  hybrid  of  computer  simulations  and  electronic 
structure  calculations  to  do  away  with  model  potentials  would  be  most  ideal  but  they  are  in  their 
infancy10  because  computationally  such  an  approach  is  so  numerically  intensive  that  presently  it  is 
restricted  to  very  short  simulated  times  and  small  number  of  atoms.  For  an  aggregate  of  even  a 
few  hundred  atoms,  the  only  feasible  study  today  is  computer  simulations  based  on  model 
potentials.  It  is  therefore  important  to  contrast  the  results  of  effective  potentials  to  these  first- 
principle  type  calculations. 
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The  overall  goal  of  the  proposed  research  is  to  develop  an  atomistic  approach  to  gain  an 
understanding  of  the  mechanisms  of  processes  to  fabricate  metal-metal/semiconductor  structures 
and  thus  to  contribute  to  their  development.  The  complexity  of  the  approach  involving  realistic 
modelling  of  atomic  interactions  and  adaption/development  of  a  suitable  dynamical  simulation 
technique  or  hybridized  techniques  can  be  broken  down  into  its  component  parts  which  can  be 
studied  stage  by  stage  and  and  can  be  grouped.  The  primary  goal  in  Phase  I  was  to  systematically 
obtain  the  energetics  and  structural  features  of  the  interfaces  of  Cu/Fe  and  to  simulate  the  dynamical 
processes  occurring  during  growth. 

Phase  I  Research  Objectives 

The  specific  objectives  of  Phase  I  follow: 

•  The  development  of  realistic  interaction  potentials  to  describe  interatomic  interactions. 
Devise  functional  forms  which  are  physically  motivated  and  can  be  parameterized  to 
experimental  data.  One  test  of  the  usefulness  of  this  approach  is  to  fit  the  potential  to  a  set 
of  experimental  /  empirical  data  and  verify  that  it  can  be  used  to  predictively  extend  these 
results  to  other  experiment  conditions. 

•  The  examination  of  energy  characteristics  of  the  heterostructure.  Conduct  energetics 
studies  of  adsorbate-substrate  structure  using  the  derived  interaction  potentials.  Examine 
strain  energy  for  lattice  mismatched  systems.  A  comparison  of  the  results  between  the 
static  and  dynamic  simulations  will  address  the  important  question  of  atomic  relaxation  and 
anharmonic  effects. 

•  The  identification  of  key  kinetic  parameters  controlling  the  growth  process  of  thin  films  of  a 
few  layers.  Simulate  the  dynamical  processes  occurring  during  growth. 

Towards  this  end.  Phase  I  laid  the  groundwork  by  delving  into  the  microscopic  energetics 
and  adatom  diffusion  studies,  dynamics  of  representative  processes  important  in  the  epitaxial  and 
superlattice  growth  phenomena.  We  proposed  to  examine  the  Fe/Cu  structure.  Although  this 
structure  is  being  studied  intensely1*4  due  to  the  novel  magnetic  properties  of  the  formed  films,  the 
emphasis  and  the  above  discussed  experimental  controversy  is  on  fee  Fe/Cu  as  opposed  to  bcc 
Fe/Cu.  Modelling  of  an  hypothetical  fee  Fe  would  not  be  realistic  because  most  of  the  bulk  data  is 
for  bcc  Fe.  It  was  also  felt  that  potential  model  -based  investigations  (that  too  for  an  hypothetical 
structure)  may  not  contribute  effectively  to  resolve  the  experimental  controversy.  Based  on  these 
observations  and  on  the  discussions  we  had  with  the  program  manager  at  the  time  of  Phase  I 
research  initiation,  we  chose  to  investigate  the  Cu/Ni  as  a  prototype  structure. 

The  rationale  for  the  choice  of  Cu/Ni  is  many  fold:  (a)  agreement  between  experimental 
investigations,1*11  (b)  lack  of  detailed  growth  simulation  studies  (the  predicted  critical  thickness 
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based  on  a  Monte  Carlo  simulation  study12  contradicts  with  the  experiment11),  (c)  availability  of 
four  different  tested  potential  models  for  the  bulk  metals,  Cu  and  Ni,  and  (d)  opportunity  to 
examine  (1)  the  effects  of  potential  models  and  (2)  the  merits  and  demerits  of  the  different 
modelling  schemes. 


RESEARCH  STATUS 


Interaction  Potential  Modelling 


Computer  simulations  based  on  atomistic  descriptions  of  the  structure  and  properties  of 
materials  require  as  input  simple  yet  realistic  interaction  potentials.  They  should  produce 
acceptably  accurate  structural  features  and  energetics  with  a  reasonable  expenditure  of  computer 
resources.  Although  most  computer  modelling  studies  to  date13  have  employed  only  pairwise 
additive  potentials,  many-body  interactions  and  directionally  nonisotropic  potentials  are  essential 
for  the  description  of  atomic  clusters,  interfaces,  and  strained  lattices.  The  two-body  potentials  are 
not  suitable  for  describing  metals/semiconductors. 

There  are  two  general  approaches  that  can  be  used  for  obtaining  potentials  for  application  in 
computer  simulation  studies.  The  first  employs  the  methods  of  quantum  chemistry  for  a  first 
principles  ( ab  initio )  calculation  on  a  small  group  of  atoms  (a  cluster)  or  local  density  functional 
(LDF)  calculation  of  the  bulk  solid.  The  second  approach  involves  fitting  parameters  in  a 
physically  motivated  function  to  relatively  limited  experimental  data  supplemented  by  results  of  ab 
initio  cluster  and/or  local  pseudopotential  calculations.  Because  of  cost,  ab  initio  methods  are 
usually  not  utilized  to  obtain  complete  potential  energy  hypersurfaces.  However,  data  for  different 
cluster  geometries  and/or  adsorbed  atom  configurations  can  be  used  to  parameterize  empirical 
potential  models  and  can  augment  often  scarce  empirical  data. 

The  embedded  atom  method  (EAM)14  was  developed  as  a  means  of  calculating  realistic 
ground-state  properties  of  metallic  systems.  We  are  drawn  to  EAM,  because  of  its  strong 
foundation  in  density  functional  theory  via  the  effective  medium  approach.  The  total  energy  is 
simply  written  as: 


E»  =  SF,(PI)  +  i  I 


IJ.  »  <J 


(D 


where  Fj(pj)  is  the  energy  needed  to  embed  atom  i  into  the  local  electron  density,  pj,  of  the 
surrounding  atoms,  and  fyj  is  a  pair  potential  between  atoms  i  and  j. 

Different  schemes15  have  been  used  to  describe  the  local  electron  density  and  the  pair 
interactions  of  Eq.  (1)  for  the  metals  Cu  and  Ni.  The  differences  between  the  various  approaches 
are  primarily  found  in  the  derivation  and  interpretation  of  the  embedding  density  of  the  electron 
surrounding  the  atom,  p ,  and  the  embedding  function  F(  p ).  Ackland  et  al15d  following  Finnis 
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and  Sinclair14*5  used  a  square-root  function  for  F(p )  and  identified  F(p )  with  the  second  moment 
of  the  density  of  states.  Whereas  the  scheme  of  Baskes  and  coworkers  scheme  (FBD),5a  is  to 
specify  p  as  an  atomic  charge  density  obtained  from  Hartree-Fock  calculations  for  the  free  atom 
and  then  to  obtain  F(  p )  by  empirical  fitting  to  the  equation  of  state  of  the  bulk.  Oh  and  Johnson's 
approach  was  to  develop  an  EAM  model  with  simple  functional  forms  which  result  in  less  number 
of  parameters  to  fit  and  thus  provide  easy  access  to  the  EAM  model  on  any  type  of  the  computer. 

The  evaluation  of  these  schemes  in  terms  of  how  the  models  perform  in  calculation  of 
unfitted  properties,  such  as  diffusion,  which  are  appropriate  for  the  planned  interfacial  and  growth 
studies  is  discussed  below  under  the  accomplishments  section  in  conjunction  with  the  kinetics. 
However,  we  feel  that  the  EAM  method  offers  a  theoretical  framework  that  holds  promise  for  the 
development  of  potentials  with  improved  transferability  from  bulk  simulations  to  crystal  growth 
and  even  cluster  chemistry.  The  trick  is  of  course  to  have  a  very  flexible  parameterization  scheme 
for  the  density,  or  embedding  term  F.  Our  studies  of  the  methodology  leads  us  to  the  belief  that 
parameterization  of  the  density  using  only  bulk  data  will  not  produce  an  EAM  method  that  is 
completely  satisfactory  for  studying  growth  processes.  However,  by  using  results  of  quantum 
mechanical  calculations  on  clusters  in  addition  to  bulk  data,  even  if  only  of  qualitative  accuracy, 
more  generally  useful  densities  can  be  obtained  that  will  satisfy  the  demands  of  simulations  of 
growth  phenomena.  This  work  must  be  deferred  until  Phase  II  due  to  time  constraints. 

Energetics  Studies 

From  a  knowledge  of  the  interaction  potential  of  the  metallic  system  under  consideration, 
the  total  energies  of  interfaces  (composed  of  the  structure  with  specific  surface  orientations  that 
differ  only  in  symmetry)  are  obtained  as  a  function  of  the  lattice  constant  ratio,  relative  orientation, 
and  relative  displacement  between  the  adsorbate  and  substrate.  These  basic  variables  can  be 
defined  in  terms  of  the  crystallographic  relationships50  between  the  two  structures  considered.  The 
substrate  model  remains  fixed.  The  variation  of  energy  as  a  function  of  the  above  variables  and 
their  values  at  the  energy  minimum  give  information  on  the  conditions  necessary  to  obtain  a 
coherent  interface. 

In  order  to  consider  the  strain  effects  due  to  lattice-mismatch  at  the  interface,  the  strain  at 
the  interface  can  be  defined  as  Strain  (%)  =  (b  -aa)/aa  ,  where  b  is  the  strained  lattice  constant 
(varying  lattice  constant  in  the  calculation)  and  aa  is  the  unstrained  (corresponding  to  the  bulk 
value)  lattice  constant  of  the  adsorbate.  The  total  energy  of  the  interface  was  calculated  as  a 
function  of  the  strain  (varying  adlayer  lattice  constant).  The  energy  curve  shows  a  minimum  with 
respect  to  strain,  which  is  dependent  on  the  number  of  atoms  in  the  adlayer. 

A  preliminary  manuscript  of  a  publication  arising  from  our  Phase  I  research  on  static 
energetics  with  simple  pairwise  additive  potentials  and  the  realistic  EAM-based  FBD  model  has 


4 


been  enclosed  as  an  appendix- 1  to  this  report  and  presents  and  discusses  our  results  in  detail.  The 
results  on  the  energetics  are  discussed  briefly  under  the  accomplishments  section  below.  We  tried 
to  address,  in  a  limited  effort  due  to  the  time  constraints,  the  important  question  of  atomic 
relaxation.  These  studies  follow. 

Atomic  Relaxation  Effects 

A  series  of  calculations  have  been  carried  out  in  order  to  examine  the  effects  of  atomic 
relaxation  at  the  dissimilar  interface  Cu  (1 1 1)/Ni(100).  Details  of  individual  runs  are  described 
below.  The  Ni  substrate  has  4  layers  with  289  atoms  in  each  layer.  The  Cu  adlayer  is  initially  a 
monolayer. 

Run  1:  Cu  adlayer  has  25  atoms  and  all  the  atoms  are  allowed  to  relax.  The  25  first  layer  Ni 
atoms  closest  to  the  Cu  atoms  are  allowed  to  relax.  All  other  substrate  atoms  were  held 
fixed  at  bulk  equilibrium  positions. 

Run  2:  Cu  adlayer  has  25  atoms  and  all  the  atoms  are  allowed  to  relax.  The  121  Ni  atoms  each 
of  the  first  and  second  layers  of  the  substrate  are  allowed  to  relax. 

Run  3:  81  Cu  atoms  in  adlayer,  all  are  allowed  to  move.  The  242  Ni  atoms  in  the  substrate  are 

allowed  to  move  (  121  in  each  of  the  first  and  second  layers  of  substrate)  as  in  run  2. 
The  initial  adlayer  is  positioned  directly  over  the  second  layer  Ni  atoms  at  a  distance  of 
1.805  A  above  the  surface  layer,  the  lattice  constant  of  the  adlayer  is  initially  3.61  A. 

Run  4:  Same  crystal  model  as  Run  3  except  initial  adlayer  has  lattice  constant  of  3.52  A  (rather 
than  3.61  A)  and  is  located  above  the  surface  at  a  distance  of  1.76  A  (rather  than  of  1.805 

A). 

Results:  It  is  clear  from  a  comparison  of  the  values  of  the  bond  length  and  the  interplanar 
distances  reported  in  table  1  that  in  Runs  1  and  2,  the  Cu  adlayer  becomes  pseudomorphic 
to  Ni.  In  Run  3,  the  initial  adlayer  "breaks  up"  into  small  clusters,  each  little  cluster, 
consisting  of  approximately  20  atoms,  is  pseudomorphic  to  Ni  substrate.  Whereas  in 
Run  4,  some  of  the  atoms  in  the  Cu  monolayer  migrated  to  a  lattice  position  above  the 
monolayer,  the  end  result  is  a  pseudomorphic  bi-layer  of  adatoms.  These  preliminary 
results  that  small  clusters  of  an  adlayer  (~  25  atoms)  will  readily  form  a  defect-free 
monolayer  on  the  substrate.  Large  clusters,  however,  do  not;  the  "break  up"  into  small 
clusters  or  form  new  layers,  and  possibly  this  is  to  relieve  the  strain  of  the  larger  cluster. 


Material: 

Cu 

Ni 

Cu  adlayer* 
Run  1  Run2 

Bond  length  (A): 

2.556 

2.489 

2.487 

2.490 

Interplanar  Distance  (A): 

1.805 

1.76 

1.754 

1.744 
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*  In  Run  I  for  the  Ni  substrate,  <  Z12  >  =  Adj2  =  -0.00778  A.  Whereas  in  Run  2,  Ni  first 
layer  Adi2  =  -0.000043  and  for  Ni  second  layer,  Ad23  =  -0.001  A 

Adatom  Diffusion  Studies 

A  knowledge  of  the  mobility  of  an  atom  on  the  surface  would  aid  in  for  assessing  and/or 
interpreting  the  order  (full  vs  partial)  of  epitaxial  growth  for  different  surface  orientations.  The 
adatom  diffusion  studies  are  also  useful  for  evaluation  of  the  potential  model  provided  reliable 
experimental  data  is  available  for  comparison  purposes  with  theoretically  predicted  values.  Instead 
of  adapting  the  cost-intensive  molecular  dynamical  simulations  for  estimating  the  diffusion 
constants,  we  used  the  variational  transition  state  theory  approach  to  compute  Arrhenius  parameters 
and  diffusion  coefficients.  Appendix-2  to  this  report,  as  a  preliminary  manuscript  of  a  publication 
arising  from  these  studies,  describes  and  discusses  our  results  on  the  self-diffusion  of  Ni  for  four 
different  EAM-based  model  potentials.  The  observed  sensitivity  of  the  variation  of  diffusion  to  the 
interaction  potential  parameters  and  comparison  with  experiment  are  discussed  below. 

Accomplishments 

Potential  Model  Evaluations  and  Self-Diffusion 

Analysis  of  various  embedded-atom  type  modelling  schemes  for  adatom  diffusion  studies, 
to  our  knowledge,  represents  the  first  study  based  on  the  realistic  bulk  models.  For  an  unbiased 
evaluation  of  EAM  models  used,  a  stringent  comparison  with  experiment  is  difficult  due  to  (a)  the 
lack  of  a  quantitative  relationship  between  the  intrinsic  (corresponding  to  an  ideal  flat  surface)  and 
mass  transfer  (involving  the  complex  processes  of  the  formation  of  surface  defects  and  their 
migration)  surface  diffusion  data  as  a  function  of  temperature,  the  former  being  consistent  with  the 
calculations,  and  (b)  the  existence  of  contradictory  data  for  the  difference  in  the  measured  intrinsic 
and  mass  transfer  data  as  well  as  large  scatter  in  the  experimental  data.  We  find  that 

•  All  the  models  exhibit  the  directional  anisotropy  for  the  self  diffusion  on  Ni  (110)  surface. 

•  At  high  temperatures,  the  models  agree  with  each  other  but  as  temperature  decreases  the 
predicted  diffusion  data  fall  within  a  scatterband. 

•  Comparison  of  the  classical  diffusion  barrier  for  Ni/Ni  (100),  Ni/Ni  (1 10),  and  Ni/Ni(l  1 1) 
indicates  that  diffusion  will  occur  more  readily  on  (111)  surface  orientation.  This  may 
imply  that  the  homoepitaxial  growth  will  be  relatively  complete  for  (111)  substrate 
orientation.  This  is  borne  out  in  a  recent  molecular  dynamics  simulation16  of  homoepitaxial 
growth  using  Lennard- Jones  model. 
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These  studies  involving  the  surface  diffusivity  and  surface  structure  with  realistic 
interaction  potentials  as  well  as  the  explosive  developments  in  the  growth  techniques  should  result 
in  renewed  impetus  for  accurate  measurements  of  adatom  diffusion. 

Energetics 

Atomistic  studies  on  static  energetics  of  Ni  monolayer  growth  on  the  Cu  substrate  have 
been  performed  using  both  central  potentials  and  the  realistic  EAM-based  FBD  model.  We  find 
that 

•  The  characterization  of  the  details  of  axial  commensuration  in  terms  of  energy  minima  and 
maxima  is  found  to  depend  on  the  atomistic  model  (rigid  vs  nonrigid,  nonrigid 
approximation  corresponds  to  the  inclusion  of  adatom-adatom  interactions)  for  the  adlayer 

•  Upon  successful  determination  of  the  embedding  energy  contribution  to  the  total  adatom- 
substrate  interaction  energy,  the  energy  charactersitics  turned  out  to  be  similar  to  those 
obtained  in  the  recent  detailed  studies  with  simple  model  potentials.  This  implies  that  for 
rigid  lattices  the  energetics  are  only  determined  by  the  relative  atomic  sizes  of  the  substrate 
and  monolayer. 

•  The  strain  energy  for  the  positive  misfit  interface  (Cu/Ni)  is  larger  than  that  of  the  negative 
misfit  interface  (Ni/Cu).  This  implies  that  a  large  number  of  layers  will  grow 
pseudomorphically  in  the  case  of  negative  misfit  interface.  This  type  of  estimate  is  of  great 
interest  in  the  fabrication  of  superlattices  made  up  of  lattice  mismatched  materials.  The 
diffusion  and  the  strain  energy  dependence  on  the  substrate  orientation  may  provide  a  clue 
to  the  choice  of  substrate  orientation  for  obtaining  better  quality  films. 

•  The  preliminary  results  based  on  the  inclusion  of  atomic  relaxations  indicate  that  the  initial 
adlayer  breaks  up  into  small  clusters,  and  each  cluster,  consisting  of  approximately  20 
atoms,  is  pseudomorphic  to  the  Ni  substrate. 
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Abstract 


Atomistic  studies  on  static  energetics  of  Cu  monolayer  growth  on  a  Ni  substrate  are 
presented.  The  sensitivity  of  the  characteristics  of  the  energetics,  as  a  function  of  the  ratio  of  the 
lattice  constants,  to  the  chosen  atomistic  model  such  as  rigid  versus  nonrigid  (nonrigid 
approximation  corresponds  to  the  inclusion  of  adatom-adatom  interactions)  for  the  adlayer  and  the 
interaction  potential  is  discussed.  Conditions  necessary  to  obtain  a  coherent  interface  between 
fcc(l  1  l)/fcc(100)  structures  are  then  examined  in  terms  of  axial  commensuration.  The  strains  and 
the  strain  energies  at  interfaces  which  exhibit  either  positive  misfit  [Cu/Ni(100)]  or  negative  misfit 
[Ni/Cu(100)]  are  obtained  employing  the  recently  developed  embedded-atom  potentials  for  Ni  and 
Cu. 
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I.  INTRODUCTION 

Extensive  experimental  work  has  been  carried  out  to  study  the  growth  of  semiconductors 
and  metals  on  a  variety  of  substrates,1  and  the  growth  of  semiconducting2  and  metallic 
superlattices.3  Theoretical  work,  until  recently,  has  mainly  focussed  on  thermodynamics  and 
kinetics  of  crystal  growth.  An  atomic  scale  understanding  of  crystal  growth  is  required  to  predict 
novel  material  properties  and  for  device  applications.  A  direct  brute-force  microscopic  approach 
would  lead  to  a  detailed  understanding  of  the  interface  formation.  Rigid  lattice  model  calculations 
of  energetics,  starting  with  the  pioneering  work  of  Frank  and  van  der  Merwe,4  have  dealt  with 
studies  of  growth  patterns,5  effect  of  substrate  induced  strains  on  the  orientation  of  monolayer 
adsorbates.6  The  advent  of  supercomputers  has  however  made  it  possible  to  undertake  detailed 
static  energetics  studies.7 

The  recent  detailed  atomistic  static  energetic  studies713  have  explored  the  dependence  of 
interface  formation  on  lattice-constant  ratio,  relative  orientation,  and  relative  displacement  between 
the  substrate  crystal  and  the  formed  Him  and  the  conditions  for  interface  coherency  of  two 
dissimilar  structures  such  as  fee  (lll)/bcc  (HO),  fee  (lll)/fcc  (100),  and  fee  ( lll)/fcc  (110). 
These  studies  are  based  on  rigid  lattice  approximation  with  pairwise  additive  potentials.  That  is  to 
say  only  the  adsorbate-substrate  (a-s)  interactions  have  been  considered  in  calculating  energy 
versus  the  lattice  mismatch  parameter,  P  defined  as  aa/as.  The  values  of  P  pertaining  to  the  minima 
and  maxima  in  the  configurational  energy,  Ua.s  has  been  interpreted  in  terms  of  axial 
commensuration  for  these  dissimilar  structures. 

In  the  calculation  of  energy  versus  P,  the  substrate  model  remains  fixed  during  the 
calculation  (in  other  words,  the  lattice  constant  of  the  substrate  does  not  change);  therefore,  the 
geometries  of  the  substrate  atoms  do  not  change.  However,  the  monolayer  adsorbate  is  changed 
during  the  calculation;  as  P  is  increased,  the  lattice  constant  of  the  adlayer  increases  and  thus 
geometries  of  the  adsorbate  atoms  change.  Thus,  the  total  energy  of  the  interface  corresponds  to  a 
sum  of  Ua.s  and  Ua„a  where  Ua.a  is  configuration^1  energ)  of  the  adsorbate-adsorbate  interaction. 
Therefore,  while  there  may  exist  a  P  which  shows  energy  minima  and  maxima  charactersitics,  the 
role  of  including  Ua.a  interactions  and  the  effects  of  atomic  relaxations  have  not  been  examined 
until  now.  The  results  of  the  latter  studies  will  be  discussed  in  a  forthcoming  publication.  Under 
such  rc  'istic  situations,  it  would  be  interesting  to  see  if,  in  fact,  the  most  stable  configuration 
occurs  when  axial  commensuration  is  attained. 

As  realistic  interaction  potentials  are  being  produced  for  metals8  and  semiconductors,9  it  is 
timely  that  quantitative  information  concerning  growth  and  interface  be  obtained  by  adapting  these 
potentials  in  computer  simulations.  Further,  it  is  hoped  that  this  quantitative  information  on 
energetics  can  be  compared  to  first-principles-total-energy  calculations  which  are  increasingly 
becoming  available.10  We  have  undertaken  such  a  study  with  realistic  potentials,  based  on  the 
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embedded  atom  method  (EAM),  for  Cu  and  Ni  metallic  structures  and  the  static  energetics  are 
considered  in  this  paper.  The  preliminary  results  based  on  the  inclusion  of  atomic  relaxations 
reported  elsewhere11  indicate  that  the  initial  adlayer  breaks  up  into  small  clusters,  each  cluster, 
consisting  of  approximately  20  atoms,  is  pseudomorphic  to  the  Ni  substrate. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2  provides  the  computational 
details  which  includes  the  lattice  models  of  the  substrate  and  adlayer  structures  used  in  the 
calculations  and  description  of  the  EAM  potential  used  in  this  study.  The  results  of  the  energetics 
and  of  the  strain  and  strain  energy  for  lattice  mismatched  interfaces  are  reported  and  discussed  in 
section  3.  Concluding  remarks  are  given  in  Section  4. 

II.  COMPUTATIONAL  DETAILS 

Lattice  Models: 

The  lattice  models  for  the  calculations  using  simple  model  pairwise  additive  potentials  are 
the  same  as  the  one  used  in  previous  energetics  studies711  for  obtaining  the  conditions  necessary  to 
obtain  a  coherent  interface  between  two  dissimilar  structures.  The  hexagonal-shaped  adlayer  was 
constructed  in  such  a  way  that  it  represents  the  first  island  formation  before  interactions  with  other 
islands  take  place.  The  substrate  consisted  of  441  atoms  per  monolayer  with  a  total  3  monolayers 
stacked  normal  to  the  growth  axis  and  the  top  layer  forms  a  two-dimensional  fcc(100)  face  with  a 
rhombic  unit  cell.  The  fee  (111)  type  adlayer  constitutes  121  atoms  on  rhombic  mesh.  The 
coordinates  of  the  central  atom  of  the  hexagonal-shaped  adlayer  is  placed  directly  on  the  top  of  the 
center  of  the  substrate  surface.  As  one  is  interested  in  investigating  the  differences  of  the  results 
between  simple  model  potentials  and  realistic  EAM  potential  models  models,  we  have  not 
considered  other  three  relative  displacements  studied  earlier  (see  fig.l  of  ref. 7b).  In  the  case  of  the 
calculations  with  EAM  model,  the  substrate  contained  289  atoms  per  layer  with  three  layers 
altogether  and  the  adlayer  is  of  121  atoms. 

For  obtaining  the  strain  and  strain  energetics  at  Cu/Ni  and  Ni/Cu  interfaces,  the  substrate 
consisting  of  289  atoms  per  layer  with  three  such  layers  and  the  (100)  orientation  and  the  adlayer 
with  (100)  orientation  has  been  constructed.  The  data  has  been  obtained  for  four  different  values 
of  the  adsorbate  atoms:  Na=25, 49,  81,  and  121. 

Interaction  Potential  Models: 

In  the  embedded  atom  method  the  configurational  energy  of  a  collection  of  atoms  is 
expressed  as  a  sum  of  a  density-dependent  embedding  energy  and  a  pair  interaction  term: 
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where  Fj(pj)  is  the  energy  needed  to  embed  atom  i  into  the  local  electron  density,  pj,  of  the 
surrounding  atoms,  and  <(>ij  is  a  pair  potential  between  atoms  i  and  j.  Different  schemes  have  been 
used  to  describe  the  local  electron  density  and  the  pair  interactions  of  Eq.  (2). 

The  EAM  was  developed  as  a  means  of  calculating  realistic  ground-state  properties  of 
metallic  systems.  Its  success  in  predicting  a  wide  variety  of  phenomena  in  metallic  systems, 
coupled  with  its  low  computational  demands,  makes  it  a  highly  desirable  alternative  to  the 
commonly  used  and  often  inadequate  pairwise  potentials  in  describing  metallic  systems.  In  the 
approach  of  Foiles,  Daw  and  Baskes  (FBD)12  which  is  adopted  in  the  present  studies,  the  local 
electron  density  pj  into  which  atom  i  is  embedded  is  approximated  by  die  sum  over  atomic  electron 

densities  pj (ry)  from  all  other  atoms  in  the  system: 


(2) 


This  functional  form  has  the  advantage  of  including  many-body  interactions  due  to  the  non-linear 
nature  of  the  embedding  function  even  though  the  computational  effort  of  evaluating  the  embedding 
energy  is  similar  to  that  for  the  pairwise-additive  interactions. 

The  functional  form  of  the  pairwise  additive  term  <J>jj(r)  is: 


4>(ry)  =  [Z(rj)Z(rj)  -Z(rcut)Z(rcut)]/r, 
=  0 


r<rcut 

r^rcut 


(3) 


where  Z^r)  is  the  effective  charge  of  atom  m  and  rcut  is  chosen  to  be  large  enough  that  the  pair 
interactions  are  very  small  at  that  value.  The  effective  charge  of  an  atom  is  given  the  functional 
form  as  follows: 


Z(r)=Z_(l  +  Prv)exp(-ar) 


(4) 


The  atomic  densities  are  approximated  by  the  single-determinant  Hartree-Fock  calculations 
of  Clementi  and  Roetti.13  For  nickel,  the  atomic  density  is  further  approximated  as  the  sum  of 
electron  densities  from  the  3d  and  4s  shells,  weighted  by  the  number  of  electrons  in  each  shell. 
We  used  the  FBD  parameters  for  all  the  calculations  reported  here. 

The  total  configurational  energy  of  the  substrate  and  the  adlayer,  within  the  EAM  model,  is 
given  by 


(5) 


where  Na  and  Ns  are  the  number  of  adlayer  and  substrate  atoms,  F  (p  )  is  the  embedding  energy 
of  the  metal  atom  under  consideration,  and  <p(r  )  is  the  short-range  pairwise  additive  interaction 
terms.  The  electron  densities  p 's  are  given  by 


N. 


p,=p; +£  p,(r„) 


(6) 


where  p!=  Xp.(r  , )  is  the  electron  density  of  the  atoms  in  the  isolated  adlayer  and  p ,  is 

i’*i  1  11 

the  electron  density  in  presence  of  the  substrate. 


N. 


pj=p;+xp1fr.i) 


(7) 


N. 

where  p]=Xp.(r  ,)  is  the  electron  density  of  the  atoms  in  the  isolated  substrate  and  p  j  is 
j  *j  1  jj 

the  electron  density  in  presence  of  the  adlayer. 


The  simple  pairwise  additive  potential,  adapted  in  the  benchmark  calculations,  are  given  by 

U(r)  =  4eexp{-|[(^)  -  l]}  (exp  {- §  [(-£)  -  l]}  -  l)  (8) 

The  parameter  o  is  fixed  according  to  the  scheme  followed  in  earlier  studies71’  excepting  that  for 
considering  the  interactions  between  adatoms  in  one  set  of  calculations,  we  have  set  e  =  1  for 
adatom-substrate  and  £=1/3  for  adatom-adatom  interactions. 


III.  RESULTS  and  DISCUSSION 


Energetics  of  fee  (1 1 1)  /  feed  00)  Structure 

The  energy  characteristics  with  and  without  the  Ua.a  interactions  as  a  function  of  the  lattice- 
constants  ratio,  P  =  aj/as,  with  fixed  relative  orientation  9  =  o°  are  shown  in  fig.  1  (panels  a’  and  a 
respectively)  using  the  pairwise  additive  potentials  of  equation  (8).  In  the  previous  studies715*  the 
details  of  the  direction  and  the  order  of  axial  commensuration  and  of  the  partial  commensuration  are 
not  only  discussed  at  length  but  also  have  been  identified  in  terms  of  the  maxima  and  minima  in  the 
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energy  seen  in  panel  a.  While  there  may  exist  a  P  which  minimizes  Ua.s,  considering  only  the  Ua.s 
interactions  does  not  minimize  the  real  total  energy  of  the  system  due  to  the  effect  of  Ua.a. 
Considering  both  Ua.s  and  Ua.a  leads  to  the  disappearance  of  energy  maxima  (at  lower  values  of  P) 
corresponding  to  the  lower  order  axial  commensuradon  in  the  x -direction. 

The  total  energy  of  the  adlayer  and  substrate  within  the  EAM  scheme  (adapting  the  EAM 
models  for  Cu  and  Ni)  has  converged  at  Z  =  0.60  as  which  is  comparable  to  Z  =  0.63  as  obtained 
with  simple  pairwise  additive  potential.  The  variation  of  energy  with  P,  using  the  EAM  model,  are 
shown  in  the  lower  most  panel  of  figure  1.  Excepting  for  the  deep  minimum  around  P  =  1,01  the 
remaining  energy  characteristics  appear  to  be  absent.  In  order  to  address  this  issue  of  lack  of 
structure  in  the  energy  curves  with  EAM  model  and  to  examine  the  dependence  of  the  energy 
characteristics  on  the  interaction  potential  model,  namely  a  simple  pairwise  additive  versus  realistic 
EAM  model,  we  first  examined  the  individual  energy  contributions  from  the  embedding  energy 
term  and  short  range  interactions  to  the  total  energy  of  the  Cu  (1 1 1)/Ni  (100)  structure. 

Using  equations  (5)-(7),  the  embedding  and  the  short-range  contributions  for  the  isolated 
adlayer  and  the  adlayer  as  well  as  for  the  total  structure  consisting  of  both  adlayer  and  substrate  as 
a  function  of  P  has  been  evaluated  and  shown  in  figure  2.  Upon  interfacing  the  substrate  with 
adlayer  and  vice-versa,  the  electron  density  increases  which  reduces  in  a  more  negative  embedding 
energy  compared  to  the  isolated  structure.  This  can  be  clearly  seen  by  examining  the  F  (p)  curves 
in  figure  2.  As  the  geometry  of  the  substrate  does  not  change  in  these  energetics  calculation  (rigid 
substrate  lattice),  the  embedding  energy  contributions  for  the  isolated  substrate  has  to  be  subtracted 
from  the  total  energy  values.  The  net  F  (p)  for  the  substrate  interfacing  with  the  adlayer  is  also 
shown  in  figure  2. 

The  short-range  repulsive  interactions  between  the  adatom-adatom  and  adatom-substrate  are 
also  shown  in  figure  2.  The  adatom-adatom  short-range  energy  contribution  to  the  total  energy  (U 
in  equation  5)  has  the  expected  behavior  as  a  function  of  P,  whereas  the  adatom-substrate  short- 
range  energy  exhibits  (not  shown  clearly  in  figure  2  but  can  be  seen  in  the  panel  c  of  figure  3 
described  below)  similar  energy  characteristics  as  that  of  panel  a  in  figure  1.  This  led  to  the 
exploration  of  obtaining  a  reasonable  estimate  of  the  contribution  from  the  embedding  energy 
interaction  in  order  to  estimate  the  total  energy  of  the  adlayer-substrate  interactions  by  adding  it  to 
the  short-range  repulsive  interaction  between  adatom  and  substrate  atoms.  Within  the  EAM 
scheme,  the  absolute  values  of  the  embedding  energy  contribution  to  such  an  interfacial  structure 
(mixed  interaction)  is  not  clearly  defined,  but  the  following  phenomenological  scheme  has  been 
adapted. 

Using  equations  (9)  and  (10),  the  total  interaction  energy  of  the  adatom-adatom  structure 
and  substrate-substrate  structure  as  a  function  of  P  is  obtained  using  separately  a  fixed  electron 
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density  appropriate  to  the  substrate  and  to  the  adlayer  corresponding  to  a  reasonable  interface 
geometry  within  the  P  space. 


i'>i 


where  p  f  is  the  electron  density  of  the  substrate  at  a  fixed  geometry  of  the  interface  (substrate 
and  adlayer),  i.e.  at  a  fixed  lattice  mismatch  as  seen  by  each  of  the  adlayer  atoms . 


j*>j 


f  • 

where  p  is  the  electron  density  of  the  adlayer  at  a  fixed  geometry  of  the  interface  (substrate  and 
adlayer),  i.e.  at  a  fixed  lattice  mismatch  which  is  felt  by  each  of  the  substrate  atoms. 

The  total  energy  of  the  adatom-substrate  interaction  is  then  obtained  by  taking  the  difference 
between  the  total  energy  (as  given  in  equation  5)  of  the  interface  under  consideration  and  that  of  the 
adatom-adatom  and  substrate-substrate  atomic  interactions  from  equations  (9)  and  (10).  This 
yields  an  expression  for  the  adsorbate-substrate  interaction 


For  a  series  of  fixed  geometries  of  the  interface,  the  above  scheme  has  been  adapted  in 
order  to  assess  the  accuracy  in  choosing  the  fixed  geometries  of  the  interface  and  the  results  for  the 
fixed  geometries  corresponding  to  P  =  0.8  and  1 .4  has  been  reported  in  figure  3.  Panels  a  and  b 
show  the  total  F  (p)  (continuous  line  curve)  contribution  (sum  of  the  first  two  terms  in  equation 
1 1)  to  the  total  adatom-substrate  interaction  energy  as  well  as  the  the  individual  energy  terms  due  to 
the  adlayer  (first  term  in  equation  11,  short -dashed  curve),  of  the  substrate  (second  term  in 
equation  12,  long-dashed  curve).  The  adatom-substrate  short-range  energy  contribution 
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(continuous  line  curve)  as  shown  in  panel  c  exhibits  energy  characteristics  similar  to  the  simple 
model  results  shown  in  figure  1,  panel  a.  Comparison  of  the  continuous  lines  in  the  three  panels  in 
figure  3  indicate  that  the  embedding  and  short-range  repulsive  energies  balance  out  in  terms  of  their 
characteristics,  i.e.,  maxima  in  one  is  a  minima  in  the  other,  and  the  magnitude  of  F  (p) ,  although 
depends  on  the  chosen  fixed  geometry,  is  much  less  than  that  of  short-range  repulsive  energy. 
Thus  the  energy  characteristics  of  the  total  adatom-substrate  interaction  energy  (short-dashed  and 
long-dashed  curves  corresponding  to  the  two  fixed  geometries)  are  same  as  that  of  those  obtained 
from  simple  pairwise  additive  potentials.  This  has  two  implications:  the  primary  implication  is  that 
the  recent  careful  examination7!5  of  the  conditions  in  terms  of  axial  commensuration  necessary  to 
obtain  a  coherent  interface  between  various  such  dissimilar  structures,  based  on  detailed  studies 
with  simple  pairwise  additive  potentials,  can  be  used  as  a  pointer.  The  second  implication  appear 
to  confirm  another  previous  simple  model  dependent  study73  which  has  shown  that  the  energy 
characteristics  are  independent  of  the  long-range  attractive  part  of  the  potential  indicating  that  for 
rigid  lattices  the  energetics  are  only  determined  by  the  relative  atomic  sizes  of  the  substrate  and 
monolayer. 

Strain  at  Cu(100)/Ni(100)  and  Ni(100)/Cu(100)  Interfaces 

The  strain  at  the  interface  can  be  defined  as  Strain  (%)  =  (b  -aa)/aa ,  where  b  is  the  strained 
lattice  constant  (varying  lattice  constant  in  the  calculation)  and  a„  is  the  unstrained  (corresponding 
to  the  bulk  value)  lattice  constant  of  the  adsorbate.  The  total  energy  of  the  interface  was  calculated 
as  a  function  of  the  strain  (varying  adalyer  lattice  constant).  The  energy  curve  shows  a  minimum 
with  respect  to  strain,  which  is  dependent  on  the  number  of  atoms  in  the  adlayer.  The  energy  vs 
strain  for  three  different  sizes  of  the  adlayer  are  shown  in  figure  4.  The  energy  difference  at  the 
zero  strain  and  the  energy  at  the  minimum  of  the  strain  is  taken  to  be  the  strain  energy.  In  our 
calculations  on  Ni/Cu  and  Cu/Ni  the  minimum  energy  was  at  strain  (%)  ~  ±  2.5  -  ±  2.6.  The 
strain  energy  as  a  function  of  number  of  adatoms  increases  monotonically  until  N  =  121  and  a 
sharp  increase  in  strain  energy  occurs.  This  might  indicate  that  clusters  with  81  or  more  atoms  will 
have  such  a  strain,  the  origin  of  which  is  not  clear.  The  question  of  how  large  must  a  cluster  be 
before  dislocations  form  must  be  investigated 

Recent  studies715  have  shown  that  the  strain  energy  as  a  function  of  the  square  root  of  the 
total  number  of  adsorbate  atoms  becomes  constant  at  some  value  of  N,  which  is  called  N(crit) 
implying  that  the  adsorbate  cluster  becomes  stable  to  the  addition  of  more  adatoms.  These  studies 
indicate  the  prediction  of  a  large  uniform  monolayer  coverage.  A  more  useful  and  pertinent 
quantity  is  probably  the  critical  size  of  atoms  above  which  the  monolayer  is  no  longer  stable  on  the 
surface  and  tends  to  force  dislocations  rather  than  continue  the  monolayer  growth.  Comparison  of 
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the  energetics  of  these  two  situations  (large  monolayer  vs  several  small  clusters  on  the  surface) 
would  give  a  nice  indication  of  what  the  system  wants  to  do. 

The  strains  and  the  energies  at  the  interfaces  of  Cu/Ni(100)  and  Ni/Cu(100)  corresponding 
to  the  positive  and  negative  misfit  respectively  are  shown  in  figure  5.  The  strain  energy  for  the 
positive  misfit  interface  (Cu/Ni)  is  larger  than  that  of  the  negative  misfit  interface  (Ni/Cu,  see 
fig.5).  Based  on  the  continuum  limit  for  the  strain  energy  in  terms  of  the  generalized  elastic 
constant,  Y,  the  thickness  (t)  of  the  overgrowth,  and  the  strain  (e),  Estrain  =  Yte  .the  large 
value  of  Y  for  the  positive  misfit  would  give  a  substantially  reduced  value  of  t  c ,  the  critical 
thickness  beyond  which  pseudomorphism  is  not  possible.  Furthermore,  a  large  number  of  layers 
are  expected  to  grow  for  a  smaller  strain  energy  per  monolayer.  Thus  it  implies  that  a  large  number 
of  layers  will  grow  pseudomorphically  in  the  case  of  negative  misfit  interface.  This  type  of 
estimate  is  of  great  interest  in  the  fabrication  of  superlattices  made  up  of  lattice  mismatched 
materials. 

IV,  SUMMARY  AND  CONCLUSION 

The  similarity  of  the  results  obtained  from  these  studies  using  realistic  interaction  potential 
model  for  the  metallic  structures  and  from  the  previous  detailed  studies  based  on  model  potentials 
suggest  that  for  rigid  lattices  the  energetics  are  determined  by  the  relative  atomic  sizes  of  the 
substrate  and  monolayer.  The  strain  energy  for  the  positive  misfit  interface  (Cu/Ni)  is  larger  than 
that  of  the  negative  misfit  interface  (Ni/Cu).  This  indicates  that  a  large  number  of  layers  will  grow 
pseudomorphically  in  the  case  of  negative  misfit  interface.  This  type  of  estimate  is  of  great  interest 
in  the  fabrication  of  superlattices  made  up  of  lattice  mismatched  materials.  The  diffusion  and  the 
strain  energy  dependence  on  the  substrate  orientation  may  provide  a  clue  to  the  choice  of  substrate 
orientation  for  obtaining  better  quality  films. 
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Figure  Captions 

Figure  1.  Energy  per  adsorbate  atom  vs  lattice  constant  ratio,  aa/as,  for  the  interfacial  structure, 

fee  (111)  /  fee  (100)  with  a  fixed  relative  orientation  of  0  =  0°.  Panels  (a)  and  (a') 
refer  to  the  results  obtained  from  using  a  traditional  central  potential.  The  panel  a’  and 
the  energies  from  the  EAM  model  for  Cu  (1 1 1)/Ni(100)  includes  the  adatom-adatom 
interactions. 

Figure  2.  The  embedding  energy,  F  ( p ) ,  per  adsorbate  atom  vs  lattice  constant  ratio,  aa/as,  for 
the  isolated  Ni  adlayer  (continuous  line  curve)  and  the  Ni  adlayer  (dashed  line  curve)  on 
the  Cu  substrate  and  the  variation  of  short-range  pairwise  interaction  energy,  and 
<p  (r ) ,  per  adsorbate  atom  (RHS  ordinate  scale)  of  the  adlayer  with  aa/as.  The  flat  dot- 
dashed  curves  refer  to  the  change  in  F  ( p )  for  the  substrate  interfaced  with  adlayer  and 
9  (r )  for  the  adatom-substrate  interactions. 

Figure  3.  Variation  of  the  individual  energy,  F  ( p )  and  9  (r ) ,  contributions  and  the  total 

configurational  energy,  U  (r)  per  adsorbate  atom  with  ag/as  for  the  adlayer-substrate 

&  s 

interactions.  Panels  a  and  b  refer  to  p  f  and  p  f  for  fixed  geometries  (see  text)  of  the 
interface  aa/as  corresponding  to  aa/as  =  0.8  and  1.4  respectively:  short-dashed  curve  is 
for  the  adlayer,  long-dashed  curve  is  the  substrate  embedding  energy  contributions  and 
continuous  curve  is  the  total  embedding  energy,  F  (  p ) ,  part  in  the  adlayer-substrate 
interaction  energy.  In  Panel  c,  the  continuous  curve  represents  the  short-range  pairwise 
interaction  energy,  <p(r ) ,  for  the  adatom-substrate  interactions,  and  short-  and  long- 
dashed  curves  correspond  to  the  the  sum  of  9  (r )  and  F  ( p )  for  the  two  fixed 
geometries  considered  respectively. 

Figure  4  The  total  configurational  energy  of  the  Cu  (100)/Ni  (100)  interface  as  a  function  of 
strain  and  as  a  function  of  the  number  of  adsorbate  atoms  (25, 49,  and  81).  The 
energies  corresponding  to  zero  strain  and  the  minimum  is  marked  with  symbols. 

Figure  5  The  total  configurational  energy  of  the  Cu  (100)/Ni  (100)  and  Ni  (100)/Cu  (100) 
interfaces  as  a  function  of  strain.  The  abscissa  and  ordinate  scales  for  the  Ni/Cu 
structure  are  on  the  top  and  RHS. 
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Abstract 


Several  potential  models  [S.  M.  Foiles,  M.  I.  Baskes,  and  M.  S.  Daw,  Phys.  Rev.  B  33, 
7983  (1986);  A.  F.  Voter  and  S.  P.  Chen,  Mater.  Res.  Soc.  Symp.82,  175  (1987);  D.  J.  Oh  and 
R.  A.  Johnson,  J.  Mater.  Res.  3,  471  (1988)  and  "  Embedded  Atom  Method  for  Close-Packed 
Metals  "Atomistic Modelling  of  Materials;  BevondPair  Potentials"  ed.  V.  Vitek  and  D.  J.  Srolovitz 
,  (Plenum  to  be  published)  World  Materials  Congress  Proc.  1988;  G.  J.  Ackland,  G.  Tichy,  V. 
Vitek,  and  M.  W.  Finnis,  Phil.  Mag.  56,  735  (1987)]  based  on  embedded-atom-method  (EAM) 
have  been  adopted  to  study  self-diffusion  on  Ni(100),  -(111)  and  -(110)  surfaces.  The  models 
consist  of  an  effective  pairwise  additive  (attractive  electrostatic  or  purely  repulsive)  and  a  many- 
body  cohesive  (correspondingly  repulsive  or  attractive)  interactions  parameterized  to  account  for  a 
majority  of  bulk  and  a  few  defect  properties.  Variational  transition  state  theory  calculations  using 
these  potentials  provide  Arrhenius  parameters  and  diffusion  coefficients,  as  well  as  a  comparison 
among  the  four  EAM  potentials  employed.  Additionally,  the  results  of  these  calculations  are  useful 
in  assessing  the  currently  parameterized  EAM  potentials  for  surface  and  interfacial  dynamical 
studies.  Discussion  of  the  results  is  given  in  the  light  of  available  experimental  data. 
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I.  INTRODUCTION 

Most  processes  occurring  on  a  crystal  surface  and  at  interfaces  are  expected  to  be 
influenced  by  the  kinetics  of  the  atoms.  The  mobility  of  atoms  on  a  crystal  surface  can  indeed  be 
orders  of  magnitude  larger  than  in  the  bulk  at  the  same  temperature  and  become  the  determining 
step  in  many  processes  such  as  surface  reconstruction,  the  rate  of  chemical  reaction,  and  epitaxial 
growth.  Diffusion  on  metals  has  received  both  experimental  and  theoretical  attention.1  Still,  little 
is  known  about  the  microscopic  mechanisms  with  which  atoms  move  on  a  crystal  surface.  Fdr 
instance,  neither  the  role  of  the  presence  of  steps  or  kinks,  nor  the  chemical  specificity  is  known. 
The  relation  between  the  surface  diffusivity  and  surface  structure  has  recently  been  undertaken  by 
detailed  molecular  dynamics  simulations  using  simple  model  potentials.  These  studies2  have 
shown  that  the  surface  structure  and  the  softness  of  the  repulsive  potential  play  a  dominant  role  in 
determining  the  diffusion  properties  and  that  complicated  migration  mechanisms,  such  as  exchange 
mechanism  for  adatom  migration  on  channelled  surface,  are  active  for  the  migration  of  adatoms. 
The  major  goal  of  this  study  is  to  critically  examine  the  use  of  embedded-atom-method  (EAM)3 
potentials,  which  have  been  fitted  to  bulk  metal  properties,  in  dynamics  calculations  of  surface  and 
interfacial  processes.  We  present  results  of  dynamical  calculations  of  diffusion  of  atomic  nickel  on 
several  nickel  surfaces  using  various  EAM  potentials. 

Variational  transition  state  theory  (VTST)4  has  been  used  successfully  in  calculating  rates 
of  gas  phase  reactions.  Additionally,  it  has  been  applied  to  other  gas-surface  diffusion  studies.5 
However,  as  is  the  case  with  all  theoretical  dynamical  studies,  the  calculation  is  accurate  only  to  the 
degree  that  the  potential  energy  surface  correctly  describes  the  system  under  study.  While 
obtaining  accurate  first-principles  potential-energy  surfaces  for  small  polyatomics  is  feasible,  the 
same  cannot  be  said  for  metallic  materials  because  of  the  magnitude  of  the  calculations.  Therefore, 
the  potential-energy  surfaces  used  in  calculations  of  metallic  bonding  interactions  usually  are  global 
functions  fit  to  reproduce  experimentally  measured  data  such  as  elastic  modulii,  equation  of  state, 
defect  energetics,  and  diffusion  rates,  etc.  Because  an  empirical  approach  must  be  taken  in 
determining  the  potential-energy  surfaces  to  use  in  interfacial  growth  related  simulations  at  the 
atomistic  level,  it  is  desirable  to  have  a  physically  motivated  model.  The  EAM  was  developed 
within  such  a  framework.  Prior  to  development  of  EAM  potentials,  the  most  prevalent  form  of  a 
potential  for  describing  a  solid  involved  pairwise  terms  only,  fit  to  structural  and  thermodynamic 
data  of  the  bulk  metal.6  However,  this  type  of  potential  is  inadequate  due  to  the  many-body  nature 
of  electron-core  interactions  inherent  in  a  metallic  material.  Daw  and  Baskes,  and  Finnis  and 
Sinclair  have  independently  developed3  a  function  which  describes  the  potential-energy  of  a 
metallic  solid  in  terms  of  the  "embedding  energy"  of  each  atom  in  the  crystal.  This  energy  is  a 
function  of  local  electron  density  that  the  atom  senses  due  to  nearby  atoms. 
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In  this  study  of  adatom  diffusion,  we  will  be  using  an  EAM  potential  which  is  empirically 
adjusted  to  static  bulk  properties.  In  addition  to  examining  the  details  of  dynamics  of  the  system, 
this  calculation  provides  an  opportunity  to  critically  evaluate  and  refine  the  EAM  potential,  as  well 
as  indicating  which  properties  of  the  system,  pertinent  to  the  dynamics  of  the  interface,  must  be 
incorporated  in  the  fit  to  correctly  describe  the  system.  The  remainder  of  the  paper  is  organized  as 
follows.  Section  2  provides  the  background  and  description  of  the  EAM  potentials  used  in  this 
study  and  Section  3  briefly  reviews  the  dynamical  methods  and  describes  the  computational 
procedures.  The  structural  models  of  the  Ni  lattices  used  in  the  calculations  are  described  in 
Section  4,  and  the  results  of  surface  self-diffusion  are  presented  in  Section  5.  Concluding  remarks 
are  given  in  Section  6. 


II.  POTENTIAL  ENERGY  SURFACES 

The  embedded  atom  method  (EAM)  was  developed  as  a  means  of  calculating  realistic 
ground-state  properties  of  metallic  systems.  The  EAM  is  based  on  density  functional  theory  and  its 
derivation  is  given  elsewhere.3  Previous  applications  of  the  EAM  to  metallic  systems  have  given 
good  results.  For  example,  energies  and  geometries,3*7  structure  of  liquid  metals,8  bulk  and 
surface  phonons,9  surface  relaxations,  reconstructions  and  other  surface  phases  of  transition 
metals,10  point  defect  properties  and  basic  alloy  properties38  calculated  with  EAM  are  in  good 
agreement  with  experiment.  Additionally,  applications  to  hydrogen-metallic  systems  have 
successfully  predicted  ordered  structures,  critical  temperatures,  and  the  influence  of  hydrogen  on 
dislocation  motion  and  fractures.3b*H  Its  success  in  predicting  a  wide  variety  of  phenomena  in 
metallic  systems,  coupled  with  its  low  computational  demands,  makes  it  a  highly  desirable 
alternative  to  the  commonly  used  and  often  inadequate  pairwise  potentials  in  describing  metallic 
systems. 

In  the  embedded  atom  method  the  energy  of  a  collection  of  atoms  is  expressed  as  a  sum  of 
a  density-dependent  embedding  energy  and  a  pair  interaction  term: 


(1) 


where  Fj(pj)  is  the  energy  needed  to  embed  atom  i  into  the  local  electron  density,  pj,  of  the 
surrounding  atoms,  and  <j>jj  is  a  pair  potential  between  atoms  i  and  j.  Different  schemes  have  been 
used  to  describe  the  local  electron  density  and  the  pair  interactions  of  Eq.  (1)  for  Ni;  they  are 
summarized  below. 
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In  the  Foiles,  Baskes  and  Daw  (FBD)  approach,  the  local  electron  density  pj  into  which 
atom  i  is  embedded  is  approximated  by  the  sum  over  atomic  electron  densities  p?  (ry)  from  all  other 
atoms  in  the  system: 


(2) 


This  functional  form  has  the  advantage  of  including  many-body  interactions  due  to  the  non-linear 
nature  of  the  embedding  function  even  though  the  computational  effort  of  evaluating  the  embedding 
energy  is  similar  to  that  for  the  pairwise-additive  interactions. 

The  functional  form  of  the  pairwise  additive  term  (J)jj(r)  is: 

Wy)  =-[Z(ri)Z(rj)-Z(rcut)Z(rcut)]/r,  r<rcut 

=  0  r>rcut  (3) 

where  Zm(r)  is  the  effective  charge  of  atom  m  and  rcut  is  chosen  to  be  large  enough  that  the  pah- 
interactions  are  very  small  at  that  value.  The  effective  charge  of  an  atom  is  given  the  functional 


form  as  follows: 

Z(r)  =  Z  D(1  +  prv)exp  (  -  ar) 


(4) 


The  atomic  densities  are  approximated  by  the  single-determinant  Hartree-Fock  calculations 
of  Clementi  and  Roetti.12  For  nickel,  the  atomic  density  is  further  approximated  as  the  sum  of 
electron  densities  from  the  3d  and  4s  shells,  weighted  by  the  number  of  electrons  in  each  shell. 
The  number  of  4s  electrons,  Ns,  is  used  as  an  adjustable  parameter  and  is  set  to  2  for  all  the 


calculations  reported  here. 

The  parameters  are  fitted  to  experimental  data  (equation  of  state,  elastic  modulii,  defect 
energetics,  data  pertaining  to  the  transition  metallic  alloys)  by  a  least  squares  procedure  and  those 
defining  the  pair  interactions  for  nickel  are  given.38  In  the  fitting  procedure,  all  of  the  parameters 
for  the  Ni-Ni  interactions  are  first  fitted13  to  the  bulk  lattice  constant,  sublimation  energy,  elastic 
constants  and  vacancy  formation  energy  with  the  exception  of  the  value  of  Ns.  For  a  given  set  of 

parameters  for  the  atomic  densities  and  pair  interaction  terms,  the  embedding  function  for  nickel  is 
adjusted  so  that  the  simplified  equation  of  state  of  the  metal.14  is  exactly  reproduced  for  successive 
values  of  the  density.  Thus  the  embedding  function  for  Ni  is  represented  numerically,  and  in  the 
present  application  is  fitted  to  a  cubic  spline 

The  differences  between  the  FBD7  and  other  approaches  are  primarily  found  in  the 
derivation  and  interpretation  of  the  embedding  density  of  the  electron  surrounding  the  atom,  p , 
and  the  embedding  function  F(p ).  Voter  and  Chen  (VC)15  have  basically  followed  the  FBD7 
scheme  excepting  for  adopting  a  functional  form  for  the  embedding  energy  (instead  of  the 
numerical  scheme  of  FBD)  and  a  Morse-type  function  for  the  short-range  interaction. 
Furthermore,  in  contrast  to  FBD,  VC  have  modelled  the  embedding  term  as  a  repulsive  interaction. 
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Oh  and  Johnson's16  approach  was  to  develop  an  EAM  model  with  simple  functional  forms  which 
result  in  fewer  parameters  to  fit  and  thus  provide  easy  access  to  the  EAM  model  on  any  type  of  the 
computer.  While  Ackland,  Tichy,  Vitek,  and  M.  W.  Finnis  (ATVF),17  following  Finnis  and 
Sinclair,36  used  a  square-root  function  for  F(  p )  and  identified  F(  p )  with  the  second  moment  of 
the  density  of  states  and  their  functional  forms  both  for  F(  p )  and  4>(r)  consist  of  a  series  expansion 
with  built  in  step-functions,  the  coefficients  of  which  have  been  fitted  to  the  experimental  data  as 
well  as  the  pressure-volume  data  from  first-principles  calculations. 


III.  COMPUTATIONAL  METHODS 


For  diffusion  on  a  surface,  in  which  all  binding  sites  are  equivalent,  the  diffusion 
coefficient  can  be  related  to  the  hopping  rate  between  adjacent  sites  provided  the  following 
assumptions  are  true: 

(1)  The  concentration  of  adsorbate  atoms  is  sufficiently  low  so  that  the  adsorbate  atoms 
do  not  interact  with  one  another; 

(2)  Diffusion  occurs  through  a  series  of  uncorrelated  hops  from  one  binding  site  to 
another  adjacent  one;  and 

(3)  The  hopping  from  one  binding  site  to  another  can  be  treated  as  an  individual 
unimolecular  process  with  thermal  rate  constant  k(T). 

This  random  walk  model  of  diffusion  gives  a  relation  between  the  thermal  diffusion 
coefficient  D(T)  and  the  unimolecular  rate  constant,  k(T),  as  follows:18 


D(T)=  2y  k<T) 


(5) 


where  a  is  the  hop  length,  and  y  is  the  dimensionality  of  the  system  (y=2  for  diffusion  on  a  surface 
wherein  all  binding  sites  are  equivalent,  y=l  for  surface  diffusion  along  channels  (such  as 
diffusion  along  the  (1 10]  direction  on  Ni(l  10)  or  along  the  [001]  direction  on  Ni(l  10)). 

The  rate  constant  k(T)  is  calculated  from  canonical  VTST  with  semiclassical  adiabatic 
ground  state  transmission  coefficients.  The  computational  procedures  used  here  are  similar  to 
those  used  for  gas-phase  reactions,4  and  details  of  the  modifications  needed  to  treat  diffusion  on 
surfaces  are  outlined  elsewhere.5  In  review,  the  procedure  consists  of  first  locating  the  reactant, 
product  and  saddle  point  geometries  (locations  at  which  the  gradient  vanishes)  then  determining  the 
reaction  path  by  following  the  negative  of  the  gradient  vector  in  mass  weighted  coordinates  from 
the  saddle  point  down  into  the  reactant  and  product  geometries.  The  reaction  coordinate  s  is 
defined  as  the  distance  along  the  reaction  path  from  the  saddle  point  and,  by  convention,  is 
negative  on  the  reactant  side  of  the  saddle  point.  The  generalized  transition  state  theory  rate  is 
expressed  in  terms  of  the  location  of  transition  state  along  the  reaction  coordinate 
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GT  kBT  QGT(T,  S) 

kGT(T,  s)  =  a  -f  -  exp(-  VMEp(s)/kBT) 

Q  (6) 

where  Qgt(T,s)  is  the  partition  function  for  the  bound  degree  of  freedom  at  the  generalized 
transition  state  location  at  s,  and  VMEp(s)  is  the  potential  along  the  reaction  path.  In  canonical 

variational  transition  state  theory  (CVT)  the  location  of  s  is  optimized  to  minimize  the  rate. 

The  generalized  transition  state  partition  functions  are  evaluated  quantum  mechanically 
using  the  independent  normal  mode  approximation  which  neglects  all  anharmonic  coupling 
between  the  modes.  The  energy  levels  for  the  individual  modes  are  computed  harmonically  unless 
otherwise  noted. 

Quantum  mechanical  effects  on  the  bound  vibrational  motion  are  included  by  treating  the 
partition  functions  quantally,  and  quantal  effects  on  reaction  coordinate  motion  are  treated  using 
multiplicative  semiclassical  ground-state  adiabatic  transmission  coefficients 

kCVT/GOD  =  k°(T)  kCVT(T)  (7) 


Within  the  adiabatic  approximation,  the  effective  barrier  along  the  reaction  path  is  obtained  from 

V°«=VMEP(S)  +  £GT(S-K  =  #>  (8) 

where  eGT(s,a=0)  is  the  ground-state  energy  level  for  the  bound  vibrational  modes  at  the 
generalized  transition  state  located  at  s.  The  transmission  coefficients  are  obtained  from  Boltzmann 
averages  of  the  probabilities  for  tunneling  through  the  ground-state  adiabatic  barrier 


kG(T)  = 


exp 


VA  /  k  T 
V  '  B l)_ 


kBT 


G  "E/kBT 
j,  dE  P  (E)  e  B 


(9) 


where  is  the  value  of  the  ground-state  barrier  at  its  maximum  and  E  is  the  total  energy.  For  the 
unimolecular  process  considered  here  the  reactant  and  product  species  correspond  to  local  wells  in 
the  adiabatic  potential  separated  by  the  adiabatic  barrier.  For  this  case,  eq.  (9)  is  only  an 
approximation  since  tunneling  does  not  occur  for  a  continuum  of  translational  energies  but  from 
discreet  energy  levels  in  the  bound  wells  of  the  adiabatic  potential 

kgd(T)  =exp(vA/kBT)lPG(Ev)  exp(-Ev/k0T) 

In  the  limit  of  a  large  number  of  bound  energy  levels  in  the  adiabatic  potential  eqs.  (9)  and  (10)  will 
be  equivalent.  However,  for  the  systems  studied,  tunneling  was  not  important. 


IV.  LATTICE  MODELS 

The  structural  model  used  in  Ni  diffusion  on  Ni(100)  consists  of  162  Ni  atoms  over  four 
layers;  41  atoms  are  in  each  of  the  1st  and  3rd  layers,  and  40  atoms  are  in  each  of  the  2nd  and  4th 
layers.  The  eight  atoms  surrounding  the  reactant  and  product  sites  were  allowed  to  move.  These 
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were  the  six  surface  atoms  surrounding  the  reactant  and  product  sites  and  the  two  2nd-layer  atoms 
lying  direcdy  beneath  the  reactant  and  product  sites. 

The  model  of  the  Ni(l  1 1)  crystal  consists  of  120  nickel  atoms  with  four  layers.  There  are 
27  atoms  in  the  first  and  fourth  layers,  and  33  atoms  are  in  the  second  and  third  layers.  The  four 
surface  atoms  and  three  second-layer  atoms  adjacent  to  the  reactant  and  product  sites  were  allowed 
to  relax. 

Two  different  models  were  used  in  the  Ni(l  10)  calculation,  because  there  are  two  different 
diffusional  directions;  one  along  the  [  IT  0]  direction  and  the  other  along  the  [001]  direction.  The 
model  for  calculation  of  nickel  diffusion  in  the  [110]  direction  consists  of  186  atoms  over  five 
layers.  There  are  42  atoms  in  the  first,  third  and  fifth  layers,  and  30  atoms  in  the  second  and 
fourth  layers.  In  this  model,  seven  atoms  were  allowed  to  relax;  the  four  surface  atoms  and  the 
three  second-layer  Ni  atoms  adjacent  to  the  reactant  and  product  sites.  The  model  for  calculation  of 
nickel  diffusion  in  the  [001]  direction  has  219  atoms,  with  49  atoms  in  each  of  first,  third  and  fifth 
layers  and  36  atoms  in  each  of  the  second  and  fourth  layers.  Atoms  (six  in  first  layer  and  four  in 
the  second  layer)  adjacent  to  reactant  and  product  sites  were  allowed  to  move.  The  number  of 
atoms  allowed  to  move  in  any  of  the  crystal  model  studied  was  not  varied  to  check  for  convergence 
of  a  diffusion  rate.  The  substrate  model  was  arbitrarily  chosen  and  used  in  each  calculation  for 
comparative  purposes. 

V.  RESULTS 

Nickel  binds  in  a  four-fold  site  above  the  Ni(100)  surface  plane,  directly  over  a  second- 
layer  nickel  atom.  Rates  of  nickel  diffusion  from  one  four-fold  binding  site  on  a  Ni(100)  surface 
to  an  adjacent  one  were  calculated  using  eq.  5  and  a  hop  length  of  2.4890  A  for  the  four  EAM 
potentials.  Comparison  of  the  calculated  diffusion  coefficients  obtained  using  the  various  EAM 
models  as  well  as  the  available  experimental  data  19  is  shown  in  Figure  1.  Additionally,  the 
calculated  Arrhenius  parameters  for  the  four  EAM  potentials  are  listed  in  Table  I. 

Nickel  can  diffuse  along  two  different  channels  on  the  (110)  surface;  the  [  11 0]  and  [001] 
directions.  The  calculated  diffusion  coefficients  are  shown  in  Figures  1  and  2  for  the  [110]  and 
in  figure  2  [001]  diffusion  together  with  the  available  experimental  data.20  Figure  2  exhibits 
experimentally  observed  directional  anisotropy  on  1 10  surfaces.  The  available  experimental  data  at 
low  temperatures  refer  to  the  intrinsic  surface  diffusion  which  is  expected  to  be  lower  than  that  of 
intrinsic  surface  diffusion.  The  Arrhenius  parameters  for  the  (1 10)  surface  are  also  litsed  in  Table 
I,  and  the  geometries  and  energetics,  reported  in  Table  II,  clearly  show  that  diffusion  will  occur 
much  more  readily  along  the  [110]  direction  than  along  the  [001]  direction,  due  to  the  lower 
diffusional  barrier  along  this  channel. 


7 


Rates  of  Ni  diffusion  on  Ni(l  1 1)  were  not  calculated  because  of  low  barriers  to  diffusion. 
As  shown  in  Table  II,  the  classical  diffusion  barrier  is  less  than  0.09  eV  for  all  EAM  potentials. 
When  the  barriers  are  too  low,  the  approximation  of  uncorrelated  hops  becomes  invalid,  and  other 
approaches  for  calculating  diffusion  coefficients  are  needed. 

VI.  SUMMARY  AND  CONCLUSION 

Calculations  of  Nickel  surface  self-diffusion  based  upon  variational  transition  state  theory 
were  performed  using  different  types  of  embedded  atom  method  potentials  in  an  effon  to  to  analyze 
the  suitability  of  the  currently  parameterized  potentials  in  calculating  surface  and  interfacial 
properties.  Stringent  comparison  with  experiment  is  difficult  due  to  (a)  the  lack  of  a  quantitative 
relationship  between  the  intrinsic  (corresponding  to  an  ideal  flat  surface)  and  mass  transfer 
(involving  the  complex  processes  of  the  formation  of  surface  defects  and  their  migration)  surface 
diffusion  data  as  a  function  of  temperature,  the  former  being  consistent  with  the  calculations,  and 
(b)  the  existence  of  contradictory  data  for  the  difference  in  the  measured  intrinsic  and  mass  transfer 
data  as  well  as  large  scatter  in  the  experimental  data.  All  the  models  exhibit  the  directional 
anisotropy  for  the  self  diffusion  on  Ni  (1 10)  surface.  At  high  temperatures,  the  models  agree  with 
each  other  but  as  temperature  decreases  the  predicted  diffusion  data  fall  within  a  scatterband. 

Comparison  of  the  classical  diffusion  barrier  for  Ni/Ni  (100),  Ni/Ni  (1 10),  and  Ni/Ni(l  1 1) 
indicates  that  diffusion  will  occur  more  readily  on  (1 1 1)  surface  orientation.  This  may  imply  that 
the  homoepitaxial  growth  will  be  relatively  complete  for  (111)  substrate  orientation.  These  studies 
involving  the  surface  diffiisivity  and  surface  structure  with  realistic  interaction  potentials  as  well  as 
the  explosive  developments  in  the  growth  techniques  should  result  in  renewed  impetus  for  accurate 
measurements  of  adatom  diffusion. 
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Table  I.  Arrhenius  Parameters 


Models 

FBD 

VC 

OJ 

ATVF 

Expt. 

Ni/Ni(lQQ) 

Ea(eV) 

0.64 

0.71 

0.45 

0.93 

1.54** 

D0(cm2/s) 

0.0038 

0.0040 

0.0047 

0.0074 

Ni/Ni(110)[110] 


Ea(eV) 

D0(cm2/s) 

Ni/NimOHOOll 

Ea(eV) 

D0(cm2/s) 


0.50 

0.0034 


1.30 

0.0049 


0.53 

0.0040 


1.67 

0.0082 


0.44 

0.0038 


1.01 

0.0061 


0.65 

0.0047 


0.65 

0.0091 


0.76*,  1.85** 


1.95*,  1.74** 


*  Table  II,  ref.  lc,  p  230. 

**  G.  Neumann  and  G.  M.  Neumann,  in  Diffusion  Monograph  Series,  "Surface  Diffusion  of 
Metals",  ed.  F.  H.  Wohlbier,  p  54. 


% 

Table  n.  Energetics  and  geometries 


FBD 

VC 

OJ 

ATVF 

Ni/Niriom 

Az  (A)1 

1.61 

1.55 

1.72 

1.64 

vt  (eV)2 

0.64 

0.71 

0.45 

0.94 

Ni/Ni(110)[ll0] 

Az  (A)1 

1.11 

1.06 

1.20 

1.14 

Vt  (eV)2 

0.50 

0.53 

0.44 

0.65 

Ni/NiCllOi  rOOll 

Az  (A)1 

1.11 

1.06 

1.20 

1.14 

Vt  (eV)2 

1.31 

1.68 

1.02 

1.54 

Ni/Nicnn 

Az  (A)1 

1.84 

1.79 

1.94 

1.89 

Vt  (eV)2 

0.06 

0.07 

0.05 

0.05 

1.  Distance  above  the  surface  in  the  reactant  binding  site. 

2.  Classical  diffusional  barrier. 
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Figure  Captions 


\ 

Figure  1.  Adatom  diffusion  on  (110)  [1101  as  well  as  on  (100)  surface  together  with  the  mass 
transfer  data. 

Figure  2.  Directional  anisotropy  of  diffusion  on  1 10  surfaces.  The  available  experimental  data  at 
low  temperatures  refer  to  the  intrinsic  surface  diffusion. 


12 


(a)m  (a)11! 


* 

V 


T  (K) 


1000/T  (1/K) 


FIG.  2 


